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Abstract 

A new procedure has been developed for the calculation of shell model coeffi¬ 
cients of fractional parentage. This procedure is free from the numerical diag- 
onalization, orthogonalization and even group theoretical antisymmetric states 
classihcation. 


I. INTRODUCTION 

New phenomena, obtained recently in nuclear physics (neutron - rich nuclei and neutron 
halo) cannot be described in shell- model approximation and require high - quality wave func¬ 
tion. One of the most efficient methods for the construction of many - particle wave functions 
with well - defined permutational symmetry and exact quantum numbers (J - total angular 
momentum and T - isospin) is the iterative procedure for the calculation of coefficients of frac¬ 
tional parentage (CFP) developed in classical works However, due to the high degree of 

degeneracy of shell - model states, the set of exact quantum numbers in many cases is not rich 
enough to classify different antisymmetrical states of system. 

The first method for classifying antisymmetrical states [§ involves the construction and 
diagonalization of the matrices representing the quadratic Casimir operators for the appropriate 
special unitary, orthogonal or symplectic groups. The real difficulty in such a case is finding 
the Casimir operator and knowing its eigenvalues. 
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For many applications the precise specification of the states is not required. It is merely 
sufficient that the states be orthogonal. The second method |Q uses the Redmond iteration 
technique in conjunction with numerical Schmidt orthogonalization or Grammian matrix di- 
agonalization. Unfortunately, although this procedure can generate complete sets for large 
spaces, the buildup of numerical errors in the process of imposing orthonormality within the 
overcomplete set of antisymmetrical states causes many problems. 

We shall here present a new and very simple method of CFP calculation, based on complete 
rejection of group - theoretical antisymmetrical shell model states classification and even of 
numerical diagonalization or orthogonalization. 

II. DEFINITIONS AND NOTATIONS 

Coefficients of fractional parentage are defined as the coefficients for the expansion of anti¬ 
symmetrical wave - function in terms of the complete set of the vector - coupled parent states 
with a lower degree of antisymmetry. In the simplest case (one particle CFP) this expansion 
is: 

(xi, . . . ,XN-l,XN\KaJTMjMT) = ^ (xi, . . . ,XAr_i;XAr| [RaJf^eN^NjN) JTMjMt) X 

KaJT 

en^nOn 

X (^KaJT] EpfipfjN II KaJT^ . 

Here one - particle variables are x* = rjcrjTj (set of corresponding radius - vector, spin and 
isospin variables) and quantum numbers ei^iji (oscillator energy, orbital momentum and total 
momentum). 

The corresponding set of many - particle state quantum numbers are K - configuration 
(it can consist of only one shell nucleons), a - all necessary additional quantum numbers and 
JT - exact quantum numbers. Semicolon means that the corresponding wave - function is 
fully antisymmetric with respect to interchanges among the nucleons variables, distributed in 
different sides of this sign. Notice the use of parentheses to represent angular momentum 
coupling. || ...) is the coefficient of the fractional parentage identificator. 
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The expression given above is written down in j — j coupling, but final results are not 
dependent on the coupling scheme. 

To calculate coefficients of fractional parentage it is enough to calculate N - particle anti- 
symmetrization operator matrix into the basis of right - hand functions (with a lower degree of 
antisymmetry). As usual, one assumes that all {N — 1) particle CFPs have been obtained at a 
previous stage. The antisymmetrizer 

-4 ^ 4..„ = ± ^ ipP 

PeSjv 

is normalized such that = A and could be presented as 

Ap...,Ar = , 

where 

K = i {1 - (/V - l)P;,_i,„} . 

Here the operator P/v-i,Ar simply interchanges the coordinates xjv-i and xjv- Because the above 
mentioned basis is antisymmetrical by xi,..., xat-i variables permutation, matrix A equals the 
matrix of a far simpler operator V: 

((KaJf,ejv£NjN) JTMjMt|F| (a'AJ'T', JTMjMt) 

+ (_l)iiv+A+'^'+'f+T+T'(iv - 1)[(2J + l)(2f + 1)(2J' + l)(2f' + 1)]5 

X E I ] I 

KaJt kJnJJ') 1 \TT' 

X l^KaJT-^eM^NjN || K'a'J'T'^'^ . 


(kajf] II KaJT) 
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A is the symmetrical real matrix - projection (A+ = A = A* , AA = A). Its eigenvalues are 
only zeros and units, and the spectral decomposition [^] is 

An,xn “Fnxrf^F 

\ / rxn 

The subscripts indicate the dimension of the corresponding matrix, n equals the dimension 
of the basis, r - rank of matrix A. Every column of F is eigenvector of matrix A, corresponding 
to a unit eigenvalue. The normalization condition of eigenvectors is 

fjr \ . . 

\ / rxn 

The matrix elements of F are coefficients of fractional parentage [^]. 

III. NEW RESULTS 

Condition Anxn-^nxn = Anxn meaus that every Anxn column is its eigenvector with unit 
eigenvalue. However, they are not normalized and even orthogonal, and only r of them are lin¬ 
early independent. The simplest method of CFP matrix calculation is based on the observation 
that the spectral decomposition of A is not dehned uniquely. The possibility exists of a free 
choice of orthogonal matrix G^xr dehned as 

A = FGG+F+ = FF+ , 

because 

F+F = G+F+FG = G+G = 1 . 

Orthogonal matrix has r(r —1)/2 free parameters, so we can choose it in a way which allows 
us to hx the corresponding number of F elements. The best choice is 

Fij = 0 if I < i < j < r . 

This means that the upper triangle of matrix F equals zero. In such a case we can obtain 
the solution of matrix equation A = FF+ in the following way. Starting with the well - known 
Redmond result P| for the hrst row 
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-^n ~ ^11 ) Pji — 1 

we can present the following columns of F in the form: 

i—1 

_ 4 , _ V- p2 

^ik 1 

k=l 

Fji = l^ij ~ ^2 

{ k=l ) 

for every value of i = 2,3, ...,r and the corresponding set of j = i + l,i + 2, Positive 

values of Fa are convenient, because the overall sign of CFP vector is arbitrary. Obviously, if 
F is constructed in this way, it fulhls the condition F’^F = 1. 

So to obtain the spectral decomposition of an antisymmetrization operator matrix with rank 
r, it is enough to calculate only the r linearly independent rows of this matrix. The calculation 
begins with the first matrix A row. If 

An = 0 , 

then every element of the first row and the first column equals zero, because by dehnition 

An = ^ Ajj . 
j=i 

This means that the corresponding basic function is not presented in fractional parentage 
expansion (it may be, that it has some hidden symmetry, wrong quantum numbers, etc.). Such 
a function has to be omitted from all the expansions in every case, when An = 0 appears. 
Because matrix A is symmetrical, every next row calculation starts from diagonal element 
{All). Then it is necessary to hnd the corresponding minor (A;xz) determinant value. If it 
equals zero, /~th row is linearly dependent on (/ — 1) the just calculated ones, because this 
determinant equals calculated rows Grammian: 

n 

Aij ^ ) An^Ajif . 
k=l 

Only in the case when An ^ 0 and det A^xz 7 ^ 0, is it worth finishing the /-th row calculation. 
Because A is a projection by dehnition, the necessary number (n) of basic states is included 
only in the case when the condition is fulhlled: 
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n M ■ 

i=i 

It is also useful to know the sum of matrix A diagonal elements, because 

Sp A = r . 

The expansions which have been presented allow after the calculation of every matrix ele¬ 
ment of A {Aij) to obtain the corresponding CFP {Fji). 

So the calculation of one linearly independent row of A gives us the coefficients for one 
antisymmetrical state expansion (one column of CFP matrix). Columns of F can be numbered 
by the positive integer numbers a = 1, 2,..., r. In this dehnition, a also equals the number of 
zeros present in the corresponding column plus one. 

IV. COMPUTATION RESULTS 

The general formalism described above was implemented in computer code to generate 
complete sets of CFP for the j = 1/2, 3/2, 5/2 shells in isospin formalism. Angular momen¬ 
tum coupling and antisymmetrization are orthogonal transformations, so is possible precise 
arithmetic with quantities no more sophisticated as square root of integer number. The CFP 
computational procedure required four different types of arithmetics: 

• algebraic and relation operations between fractions composed of 16 - bit integer numbers 
for operating with quantities like angular momentum. 

• square root of 64 - bit integer numbers. 

• algebraic and relation operations between numbers represented in the form n/m/^/k, 
where n,m,k are 64 - bit integer numbers for both middle scale and 6j coefficient calcula¬ 
tions. 

• algebraic and relation operations between numbers represented in the form n/m/^/k, 
where n,m,k are arbitrary length integer numbers accomplished as variable long vector 
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with prescribed position weight for large scale calculations in evaluation A and F matrices. 
Orbit j = 5/2 was calculated with 20 - position vectors and position weight equal to 1000. 

As a preliminary step complete enumeration of the states in isospin formalism was carried 
out. By means of combinatorial calculations all JT and number of states with the same JT 
were found. This enable hud out all permited CFP and to characterized them uniquely by 
means of parent and daugther states. Results for the j = 5/2 nuclear shell are presented in the 
table. 


N 

T 

n(J)m 

2 

0 

1(1)1 1(3)1 1(5)1 


1 

1(0)1 1(2)1 1(4)1 

3 

1/2 

2(1/2)1 4(3/2)l 6(5/2)2 5(7/2)2 4(9/2)l 3(11/2)1 2(13/2)1 


3/2 

2(3/2)l 3(5/2)l 2(9/2)l 

4 

0 

2(0)2 7(2)3 8(3)1 8(4)3 7(5)1 5(6)2 2(8)1 


1 

7(1)2 10(2)2 11(3)3 11(4)2 9(5)2 6(6)1 4(7)1 


2 

1(0)1 3(2)1 3(4)1 

5 

1/2 

9(l/2)2 16(3/2)3 21(5/2)4 22(7/2)4 21(9/2)4 17(11/2)3 13(13/2)2 



8(15/2)1 5(17/2)1 


3/2 

6(1/2)1 11(3/2)2 14(5/2)2 14(7/2)2 13(9/2)2 10(11/2)1 7(13/2)1 


5/2 

3(5/2)l 

6 

0 

11(1)3 17(2)1 20(3)5 20(4)2 18(5)3 15(6)2 11(7)2 



4(9)1 


1 

6(0)2 17(1)1 26(2)5 30(3)3 30(4)5 26(5)2 21(6)3 



15(7)1 9(8)1 


2 

7(1)1 10(2)1 11(3)1 11(4)1 9(5)1 


3 

1 (0)1 
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1/2 

14(1/2)2 25(3/2)3 32(5/2)4 35(7/2)4 34(9/2)4 29(11/2)3 22(13/2)2 



15(15/2)1 10(17/2)1 


3/2 

10(1/2)1 18(3/2)2 23(5/2)2 24(7/2)2 23(9/2)2 18(11/2)1 14(13/2)1 


5/2 

6(5/2)l 

8 

0 

4(0)2 17(2)3 20(3)1 20(4)3 18(5)1 15(6)2 7(8)1 


1 

17(1)2 26(2)2 30(3)3 30(4)2 26(5)2 21(6)1 15(7)1 


2 

3(0)1 10(2)1 11(4)1 

9 

1/2 

9(1/2)1 16(3/2)1 21(5/2)2 22(7/2)2 21(9/2)1 17(11/2)1 13(13/2)1 


3/2 

11(3/2)1 14(5/2)1 13(9/2)1 

10 

0 

5(1)1 8(3)1 7(5)1 


1 

3(0)1 10(2)1 11(4)1 

11 

1/2 

6(5/2)l 

12 

0 

1 (0)1 


First column contains nucleon number N in the shell, second stands for total isospin T, 
last column gives total J and two additional numbers: n - CFP number for specihed JT (the 
dimension of corresponding A matrix); m - stands for number of states with the same JT (the 
rank of a corresponding A matrix). Thus we can have n x m CFP with the same JT. Though 
generally j = 5/2 nuclear shell has 3359 CFP, computational method determines m{m — l)/2 
zeros by zeroing F matrix upper right corner. This gives 141 zeros in complete CFP set. Actual 
number of CFP equal to zero is 328 due possible nonincounted symmetries. 

To illustrate the evaluation of the CFP consider N=4 state with J = 6 , T = 0. As follows 
from the combinatorial calculation in this case A matrix dimension is 5 and rank 2. 


A [6,0] _ 

^4 — 


7 

14 

-21 

28 

12P33 

3P462 

2v^2002 

P6006 

53 

-56 

-147 

14 

396 

99v44 

22^546 

33P182 


49 

21 

35 


99 

11P39 

33vT3 



175 

-21 



286 

143 vA 


49 

143 J 
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Computed matrix is matrix projection, its rows and columns are labeled by CFP left - hand 
bracket indices. A matrix spectral decomposition gives F matrix, which has 2 columns(as 
predicted by combinatorial calculation) and contains CFP as matrix elements. To facilitate 
comprehension F matrix is displayed in one column. 

(3(2): 7/2,l/2;5/2,l/2||4(2):6,0) 

(3(1): 7/2,l/2;5/2,l/2||4(2):6,0) 

(3(1): 9/2,l/2;5/2,l/2||4(2):6,0) 
(3(l):ll/2,l/2;5/2,l/2||4(2):6,0) 

(3(l):13/2,l/2;5/2,l/2||4(2):6,0) 

(3(2): 7/2,l/2;5/2,l/2||4(l):6,0) 

(3(1): 7/2,l/2;5/2,l/2||4(l):6,0) 

(3(1): 9/2,l/2;5/2,l/2||4(l):6,0) 
(3(l):ll/2,l/2;5/2,l/2||4(l):6,0) 
(3(l):13/2,l/2;5/2,l/2||4(l):6,0) 

CFP are denoted (iV(d) : J, T; j, t||iV(Q;) : J, T). First number in the bracket is iV — 1, 
second number closed in parentheses stands for parent state number with the same total angular 
momentum J and isospin T, shell’s j and t stands after semicolon. The same is for daughter 
state in the right - hand. Columns of matrix F form orthonormalized sets of CFP: columns 
with the same a are normalized to unit and orthogonal. 

V. CONCLUSIONS 

The suggested method of CFP calculation is not dependent on a coupling scheme and con¬ 
figuration complexity. It is the same in every case. So, it is useful in the description of highly 
excited states, when we try to investigate more interesting phenomena. The method is useful 
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27i5 

7 

6^55 

28 

3^770 

-63 

TsoMo 

56 

v/TooTo 


6 

7330 

-21 

Tnss 

-49 

THo^ 

-21 

TTMIs 


9 



for intrinsic CFP calculation, when the system centrum of mass state is under control in ev¬ 
ery antisymmetrical basic function. It is for this reason that we have introduced the original 
CFP abbreviation and dehned the procedure for arbitrary conhguration. How about possible 
numerical instabilities, we have developed arithmetic, which allows calculations without approx¬ 
imations. Numbers are represented in the form: n/m/y/k, where n,m,k are integer numbers, 
and all calculations are made with numbers in this form. Our computational experience shows 
that method has high accuracy, therefore real numbers are insufficient. We have used 25 Mhz 
AT - 386 for calculation of CFP of the nuclear shell j = 5/2 and it take us about 20 hours. 

This work was supported, in part, by a Soros Foundation Grant awarded by the American 
Physical Society. 
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